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Abstract 

We obtain exact upper and lower bounds on the steady state drift velocity, and 
kinetic energy of electrons, driven by an external field in a weakly ionized plasma (swarm 
approximation). The scattering is assumed to be elastic with simplified velocity dependence 
of the collision cross sections. When the field is large the bounds are close to each other and 
to the results obtained from the conventional approximation of the Boltzmann equation 
in which one keeps only the first two terms of a Legendre expansion. The bounds prove 
rigorously that it is possible to increase the electron mobility by the addition of suitably 
chosen scatterers to the system as predicted by the Druyvesteyn approximation and found 
in experiments. 

PACS numbers: 52.25.Fi, 05.60,+w, 52.20.Fs, 02.30.Mv 

I. Introduction 

The behavior of the electron mobility in a gas composed of several species is a subject 
of continued experimental and theoretical investigations [1-4]. Of particular interest is 
the fact that the addition of certain types of scatterers, i.e. neutral species, to the gas 
increases the electron mobility and therefore the electron current in an applied electric 
field [3,4]. This effect is potentially of practical utility and, as was pointed out by Nagpal 
and Garscadden [4] , can be used to obtain information about scattering cross sections and 
level structure of different species. 

The fact that the mobility can actually increase with the addition of scatterers is at 
first surprising: it is contrary to the well known Matthiessen rule in metals which states 
that the total resistivity due to different types of scatterers is the sum of resistivities due 



1 



to each of them [5]. A closer inspection shows that Matthiessen's rule refers to the linear 
regime of small electric fields while the observations and analysis in gases [3,4] are in the 
nonlinear high field regime. 

This still leaves open the question of the validity of approximations commonly made in 
calculating the current of weakly ionized plasmas in strong fields. We therefore investigate 
here rigorously the stationary solutions of the kinetic equation for the electron velocity 
distribution function in cases where the electron-neutral (e-n) collisions are purely elastic 
and their cross section is modeled by a simple power dependence on the electron speed. 
In particular we establish two-sided bounds for the electron mean energy and drift in 
the presence of an external electric field. These bounds show that the results obtained 
for the current and energy of the electrons in the usual approximation, which neglects 
higher order terms in a Legendre polynomial expansion and gives the Druyvesteyn-like 
distribution for large fields, are qualitatively right and even provide good quantitative 
answers. In fact they are sufficiently precise to confirm an increase in the current for large 
(but not for small) fields upon addition of some gases, provided the mass of the added 
species is smaller than that of the dominant one, e.g. adding Helium to an Xenon gas, and 
the different cross sections satisfy certain conditions. We believe that our analysis can be 
extended to include more realistic elastic cross sections and inelastic collisions; these are 
most important in practice for enhancement of the electron mobility. 

II. Kinetic Equation 

Our starting point is the commonly used swarm approximation, applicable to gases 
with a very small degree of ionization, [6-10]. In this approximation only e-n collisions 
are taken into account in the kinetic equation for the electron distribution function (EDF) 
/(r, v, t). The neutrals themselves, which may consist of several species, are assumed 
to have a Maxwellian distribution with a specified common temperature T n . Further 
simplification is achieved if the e-n collisions are assumed to be essentially elastic: the 
collision integral can then be reduced [1,6] to a differential operator due to the great 
difference in the masses of the electrons and neutrals. To simplify matters further we 
consider the case where the scattering is spherically symmetric. The stationary kinetic 
equation for the normalized EDF, in a spatially uniform system with constant density n 
subject to an external electric field F, can then be written in the form [6] 



F-V v f=-— e(v)—- 

m v z ov I a{v) 



( 
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mv dv J X(v) 
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Here e,m are the electron charge and mass, <7j is the collision cross section with species i 
whose mass is Mj and number density is N i: A is the mean free path in the e-n collisions, k 
is Boltzmann's constant, fo is the spherically symmetric part of the distribution function, 



fo(v) = ^J /(v)rfO. 



We note that e is a small parameter equal to the ratio of the electron mass to the mean 
mass of neutral scatterers, e = ml" 1 , where M _1 = M~ 1 N i a i / N^i. 

A. Velocity independent cross sections 

We shall consider first the case where (Ti(v) is independent of v so A = const and 
e = const. Taking the electric field parallel to the z-axis Eq.(l) can be written in the 
following dimensionless form 



du z u 2 du 



4 / f Tdf 
u [ fo + -7T- 

u ou 



+ u(f -f), (2) 



where 



u = 7 v, u = \/u£ + u$ + u z , 



7 = t/^, T=% E pA|F| 



kT ' T ' kT 

with some fixed T specifying the units of the temperature. We normalize / so that 

i- J f(u)d 3 u = J™ u 2 f du = 1, (3) 

When E = the stationary distribution is the Maxwellian with temperature T, 

/ = /o = M(.) = yXexp(^), (4) 

M(u) is the unique solution of (2) for E = 0, e ^ 0. When E ^ the situation is more 
complicated. Only for E small compared to e can we expect the stationary EDF to be 
close to M(u). But in the physically interesting regimes it is e which is small compared 
to E. On the other hand if e ~ the collisions almost do not change the electron energy 
so it is difficult for the electrons to get rid of the energy they acquire from the field. The 
limit e — > is therefore singular. In particular there is no well defined reference stationary 
state for e = about which to expand the solution of (2). 

B. Legendre Expansion 
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The usual method [8] of solving (2) is to expand /(u) in terms of the Legendre 
polynomials Pi, 



00 21 + 1 f 

/(u) = Y J fi{u)Pi{cos9), /,(«) = / f(u)P l (cose)dQ u , 

1=0 71 J 



(5) 



where 9 is the angle between u and the field F: cosQ = u z /u. Substituting (5) into (2) we 
obtain an infinite set of coupled ordinary differential equations for / > 0, u > 0. These 
have the form 



E (dh 2 \ Id 



3 \ du u' 



1 1 , . T df 
u [fo + —- r 

u du 



I = 0, 



(6) 



and 
E 



21 



_ / dfl-! _ l ~ 1 fi _\ + 1 + 1 ( <Vl+l + 1 + 2 f 

1 \ du u ) 21 + 3 \ du u 



uf h 1 = 1,2,.... (7) 



Eq.(6) can be integrated to give, 



3e 



^ E Y^^udu 



Tdfo 



(8) 



where the arbitrary constant of integration was taken to be 0, using reasonable assumptions 
on the behavior of / as u — > and u — > 00. 

In the conventional [8-10] approximation scheme only two terms of expansion (5) are 
kept. This is equivalent to assuming fi(v) = for / > 2. One then adds to (8) one more 
differential equation, obtained from (7), for / = 1 



E— =uf ± . 
du 



(9a) 



Substituting (8) into (9) then yields an equation for f 

3e7 

^2 



, . 3eT 2 \ df 3e 3 



whose solution is 



/o = Cexp - 



x 3 dx 



(96) 



(9c) 



/ Tx 2 + £; 2 /3e 

This /o becomes the Maxwellian M(u), (4), when E = and the Druyvesteyn [11] 
distribution / D when T = 0: 



(10a) 
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where T is the gamma function. Using (9) and (10a) one can find f\\ 



/l = _C?f!ex P (-^). (106) 

For T > 0, f Q in (9c) will always have a Maxwellian form for u » (E 2 /Te) 1/2 . 

The first two harmonics are sufficient to find the mean energy per particle W and 
mean speed (drift) w of the electrons which are physically the most important properties 
of the stationary state, 

W=-J v*f(v)d\ = w J o u*f„du, n, = — 2 J uji'u = _ Jf ,/ h du. (11) 

We shall now study the properties of these moments without the approximations made for 
explicitly solving Eq.(2). 

III. Moments of the Distribution Function 

We assume that moments 



(0 



/>oo 

= / u^^du (12) 
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exist at least for < k < 9. Multiplying (7) by a positive power k of u and integrating 
over u yields the equations 



E 



J + k-l (,_!) (I + !)(< + 2 -AQ m +1) 
1 2/-1 • M& - 1 + 2TT3 • Mfc - 1 



^2.1, (13) 



In terms of these moments w and can be written, using (11) and (8), as 



e 

w = — — 

E'y 



Mf ] - ATMf ] 



, W=^M?l (14) 



We will now construct estimates of w and W by using (8) and (13) to get relations between 
the M^. 

i) Taking / = 1 and k = 3 in (13) and substituting (8) for the calculation of gives 

M 0, = i = ^(A<r-5rx« 0) ). (15) 

ii) For / = 1, k = 6, Eqs. (13) and (8) yield 

M^ + ImW = ^-(mP-STM?\ (16) 
5 



iii) The set / = 2, k = 4 allows us to find M 5 
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(2). 



Mf ] = -i-EM^ = 10e(M 0) - 4TA<< >) 
and eliminate it from (16) to obtain, 

(1 + 2e)A4 0) " ^TeMf = ^(M ( 9 0) - STM^). (17) 

Further calculation using different / and k will give additional equations for the M^ which 
might improve the estimates, but we shall use here only (15) and (17). 

Exploiting now general bounds on moments of the nonnegative density fo(u) derived 
in the Appendix we obtain two-sided bounds for M^\ M 5 , which determine, by 

(14), the electron drift w and mean energy W. 

Inequalities 

The upper bounds on Mj, j = 3,4, 5, (we have dropped the superscript zero) can be 
calculated from (15) using (A5): 



M 4 < M\ /z => 1 > (Ml - 5TM 4 ) => Ml - 5TM± < 0. 

By solving the last inequality one gets 



5T IE 2 /5T X 2 



M, < a. a — ■ y — ■ [~ ) ■ (18) 
The same technique using bounds, 

-M 3 <(^ 6 ) 1/4 , M 5 <(M 6 ) 3 / 4 

gives 

M s < a 1/2 , M 5 < a 3/2 , M 6 < a 2 , ^ > a. (19) 

M4 

The derivation of lower bounds via (15) and (17) is more intricate. Keeping in mind 
that e is small, we use (17) in the form of an inequality 



^ (1 + 2e)> A^" 8T A^-V^ V^" 8T| ' 



where we have used M7 < ^M^Mg in virtue of (A5). Using now (A6) with j = 5, n 
1, s = 4 we obtain 



M 9 . ( M 6 S 



M 5 ~ \M 



15 



6 



and a quadratic inequality for Mq/ M. 5 whose solution is 

£s<*/», 6 = 4r + ,/(4r) 2 + ^ 2(1 + 2t) (20). 

.M5 V e 

We repeat now in (20) the use of (A6) with i = 6, fc = 1, s = 2 and i = 6, k = 2, s = 3/2 
with the results 

^< 6 , (21) 
M 4 M 3 

One can solve (15) for Me in terms of M4 and using (21) obtain the inequality 
M 4 = = ^ f ^ + 5TM 4 ) > 6" 1 f ^ + sW) . 

Its solution is 

m * a it^W- p2) 

Similarly expressing .M5 and .M3 through Ai^/Aie and .M3/.M6 respectively and using 
(15), (20)- (22) we find the lower bounds. Together with (19) they allow us to write down 
two-sided bounds for Aij, (j = 3, 4, 5) in the form 

am ., > Mj > (23) 

These are sufficient, by (14), for the estimation of w and W. One can write immediately 

ma m-E 2 . 

W- W - W-5D - (24a) 

Using the definition (14) and the inequality (A5) we obtain 

> w > ±.My\M\ lz - 4T), (246) 

which can be combined with (23) for j = 5 to get explicit bounds on w. 

The lower bounds in (23) are useless when E — > and the solution of (2) approaches 
the Maxwellian. Generally, the inequalities (23) become more useful the larger E is. 

IV. Comparison with the Druyvesteyn Approximation 

When the background temperature T is small compared with i?e -1 / 2 it can be 
neglected in (18), (20) and the bounds (24) look simpler: 



e 



/ 4 VE e 1 / 4 ^ mE mE 



>™> -7KT, , o-mim . T^=>^> - 7= ?= , = F - (25) 



7 " - 7[2(l + 2e)]V4' 2 7 2 ^- " 2 1 ^2e(l + 2e) 
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These bounds specify the electron drift and mean energy as functions of the electric field 
and gas parameters within errors of about ±20% for the mean energy and ±8% for the 
drift uniformly in E and e. For comparison w and W obtained from the Druyvesteyn 
distribution [10a] are 



w « 0.897- 



e i/4£i/2 mE 



W « 0.854 



7 



27 2 v^ 



(26) 



in good agreement with (25) when e << 1. 

Experimentalists also measure sometimes the transversal D t and longitudinal Dl 
diffusion constants for the electron swarm. While Dl cannot generally be expressed [2,9] 
in terms of the velocity moments, 

D t = D= ^M 3 

37 

is just the isotropic diffusion constant, where A is the mean free path of electrons (A = A 
here). When T can be neglected we obtain 



A 

37 



1/4 



>D> [2(l + 2e)]- 3 / 4 A 



1/4 



(27). 



For comparison 



A (E*\ 1/4 
D « 0.759— — 
37 V e 



in the Druyvesteyn approximation. 

V. Mobility in Binary Mixtures 

The increase of electron mobility w/F in a plasma upon the addition of a small 
amount of a new gas has been observed in [3] . It was calculated in [4] within the two-term 
approximation (8), (9) for binary mixtures of a heavy noble Ramsauer gas and Helium 
addition. We shall show here rigorously that this effect exists even with constant collision 
cross sections. Using (11) gives 

(28) 



37 3 



and for / = 1 Eq.(13) reads 



M 



(i) 
fc+i 



E ( -kMfl, + 2 



fc-i 



(29) 
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When E — > we may neglect the second term in (29) and obtain 

2E . (o) 4E 2 [2 eFX 



w ~ — m == = -\ ; , (30) 

using (4) and the initial notation. The resistivity F/enw is here proportional to Yl Ni<Ji, 
which is just Matthiessen's rule. 

Let us consider now the case of a strong field, kT << eFX/ ^Je, for a binary mixture 
i = 1, 2 and use the two-term ansatz (8), (9). We then have the Druyvesteyn distribution 
(10) with the moments (26). Using (14) and the notation 

N 2 M X (T 2 

a = , u = , 9 = — 

iVi + iV M 2 ' en 

we can write explicit expressions for the drift and mean electron energy 



w = 0.897 J- ^ 31 

l/ (Nt +N 2 )a 1 ^^M[ (1-a + aO 3 / 4 ' v ; 



W = 0.427 — J - a + a9)(l-a + a^)]" 1 / 4 . (32) 

(Wi + 7V2J V m 

Both the current and energy of electrons increase, but the mobility w/F decreases, as the 
field F increases.. 

Let us now keep the total gas density N\ + N 2 constant and vary the relative 
concentration of components by changing a. A simple analysis of (31) shows that w 
can be no n- monotone when both 9 and \x are larger than 1. For example, if 9 = 5, \i = 20 
then considering w as a function of a, w = w(a), we have 

%^1.41, 44 ^ 0.95 
w(0) w(0) 

Here w(a m ) is the maximum value of w obtained for a m ~ 0.11. The drift speed is almost 
the same in the pure species 1 and 2, but it is noticeably larger in a mixture. The mean 
energy of electrons changes more: when the lighter component substitutes for the heavier 
one it goes down: 

W(0) W(0) 

There is even a more striking situation, when one just adds the lighter gas keeping 
the density Ni of the heavier component constant. In this case 

W{6) ~ % + +60)V^ W{6) ~ (1 + 6) ~ 1/2{1 + 69) ~ 1/4{1 + ^ )_V4 ' (33) 
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where 5 = N 2 /Ni. Increasing 5 we increase the density of scatterers, but for 5 = 5 m = 8.5% 



W(5m) 

w(0) 



1.4, 



while the electron energy decreases: W(5 rn ) w 0.5W(0). 

We obtained these results approximately - by truncating the series (5). However 
comparing (26) with the bounds (24) we see that the drift velocity and mean energy for 
the Druyvesteyn approximation cannot differ from the exact solution by more than about 
+12, —6% and ±17% respectively. Hence the non-monotone dependence of the electron 
mobility on the density of the light species holds for the exact solution of the kinetic 
equation (2). When we had w max ~ 1.40mj(0) (within the approximation) a possible 
exaggeration of w max by 12% and underestimation of w(0) at most in 6% could reduce 
their ratio from 1.40 to 1.16 but the effect is clearly there without approximations. 

The explanation of such unusual behavior of the electron drift in the nonlinear regime 
is quite simple. When M 2 < Mi the addition of species 2 makes the energy transfer from 
the electrons to atoms easier in the elastic collisions. Consequently the mean electron 
energy W will drop leading to a net increase of the mean free time r(v) ~ X/v. The 
competition of A and v is shown by formulas (31) and (33) where a, 5 represent the 
concentration of the lighter species and \i is proportional to its relative effectiveness in 
the energy transfer. Adding about 10% of a component with atoms of mass rri2 ~ 0.05mi 
the mean electron energy decreases by about 1/2 implying the increase of w by about 40%. 

This rise of the electron mobility can be stronger [4] in the case when the collision 
cross section of the main (heavy) component is energy dependent and decreases with the 
electron energy. 

VI. Simple velocity dependent collision cross sections 

We consider here a one-species plasma with the atoms of mass M and generalize the 
bounds (24) for e-n collision cross section of the form 




(34) 



where the exponent p can be positive or negative in a certain range. Setting 




eF 



kT, 



n 



rn 



mNo~o ' 



mv, 




we can rewrite (1) as 



-e 4 



p+2 df 



dy z 



y 2 dy 



1 d 



y 




y dy ) _ 



+ y p+1 (fo-f), 



(35) 
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where v = e~ 1 / 4 v y and we have in mind situations with "strong" electric field t « 1. 
Using the Legendre series expansion (5) for /(y) we again obtain the infinite set of coupled 
equations for harmonics fi(y) 



4 y +p fi = 



i 



21 + 1 V dy 



dfi 



i-i 



l-l 



y 



fi-i + 



i + i fdf i+1 1 + 2 



21 + 3 \ dy 



+ 



y 



-fi+i 



(36) 



for I = 1,2, 3, ... and one more equation 



fi = -3e 



fo + 



t_dfo 
y dy 



(37) 



corresponding to (8). 

Methods similar to those in the Section 2 allow us to derive the pair of equations for 
moments, which generalize (16) and (17): 



where 



M(2p + 6) = e p/2 M(2), M(3p + 9) = e p/2 cM(p + 5), 

1 f°° 
c=-[p + 6 + 4e(p + 3)], Mk) = j o fo(y)y k dy 



(38) 



and the background temperature parameter t is neglected for simplicity. In terms of these 
moments, which clearly satisfy (A2), we have for the electron drift and mean energy 



i-p 



w = e 4 v M(p + 5), W = e 



(39) 



A calculation similar to that described in Section 2 and Appendix shows that 
Eqs(38),(39) yield the following upper (U) and lower (L) bounds for for w and W: 



wl <w < wu, Wl < W < Wjj, 



W L = Vq 



p+1 
e \ 2 P +4 



p+i 



Wjj = Vq€ 2 p+ 4 , Wl 



mvr 



-e p+ 2 c p+ 2 



W, 



mvr 



u 



-e p+ 2 



2 2 
which give (24) for the velocity independent cross section p = when T « e _1 / 2 i?. 
We can find the approximate solution of (35) 



(40) 



/ D (y) = Cexp 



x 



2p+3 d 



X 



o e p/ 2 + 3tx 2 + 2 P 



(41) 



using the two-term ansatz which leads to the Druyvesteyn function (9c) for p = 0. 
Computing the moments in (39) with the help of (41) yields the explicit formulas 



p+i 



wd = € 2 p+ 4 vq 



P+3 



~2p + A~ 


2p+4 / 


3 





_p+6 
> + 4 



2p + 4 ; ' 
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The bounds in (40) for the drift and energy as functions of the parameter p are shown 
in Fig.l in the form wb/wd — 1, Wb/Wd — 1 respectively with the Druyvesteyn result 
(42) for comparison (we use the subscript "B" for both "L" and "U"). The accuracy of 
two-term approximation for our models is quite good. 



Wd = e p+ 2 



mvr 



2p + 4 
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Appendix 

The moments Aik involved in (15), (17)-(24) are the integrals of the non-negative 
function fo(u): 

foiu) = lj f^ sm6de - 

We can easily show that \rvM.{k) is a concave function if one treats fcasa continuous 
variable: 

£*M>0. (Al) 

(Al) is equivalent to the inequality 
which can be written using (12) as 

oo poo / poo \ 2 

x k fo(x)dx-j y k In 2 (y)f (y)dy- I J x k \nxf (x)dx) = 

\ 1^ 1^ x k y k In 2 f (x)fo(y)dxdy > 0. 
The concavity implies obviously 

k — % n — m 

Taking k — i = n — m, n — k = j we obtain 
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For the case k = m (A3) yields inequality 

(M k y~ % < (M i ) j - k (M j ) k -\ 0<i<k<j, (A5) 

which is equivalent to the following useful set: 

' M j+n \ < M 

7+sn 

/ Mi \ M 

. Mj ) ~ 
where i,j,n, k>0, s > 1 and % > sk. 



Mj J ~ Mj ' V-Mi-J ~~ Mi- S k ' 1 J 
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Figure caption 

Fig.l. The bounds of the electron drift (Fig. la) and mean energy (Fig. lb) as functions 
of exponent p in (34) . 
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